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NOMENCLATURE 

harmonic well optical potential parameter (central te’.’m), 
dimensionless 

harmonic well optical potential parameter (central term), fm*^ 

harmonic well optical potential parameter (spin-orbit term), 
dimensionless 

harmonic well optical potential parameter (spin-orbit term), fm“^ 
impact parameter vector, fm 

projectile kinetic energy in center of mass system, MeV 

proton charge, esu 

defined in equation (7) 

spin non-flip scattering amplitude, fm 

spin-flip scattering amplitude, fm 

Planck's constant, jniile-socond 

total angular momentum v^octnr, dimensionless 

center of mass momentum wave number, fm"^ 

orbital angular momentum voctor, dimensionless 

nucleon mass, kg 

unit vector normal to scattering plane, dimensionless 

Legendre polynomial of order L, dimensionless 

associated Legendre function of first kind and order L, 
dimensionless 

incident beam polarisation vector, dimensionless 

relative position vector, fm 

WKB turning point, fm 

coulomb turning point, fm 

spin angular momentum vector, dimensionless 

ii 
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U* defined in equation (22) 

Uc defined in equation (21) 

Vc central potential term, MeV 

Vs spin-orbit potential term, MeV 

+ 

Veff defined in equation (23) 

V relative velocity in center of mass system, m/sec 

Zj target nuclear charge number, dimensionless 

B defined in equation (11) 

r gamma function of complex argument, dimensionless 

6l“ defined in equations (16) and (20) 

n|_ defined in equation (13) 

a twice the spin angular momentum, dimensionless 

0 ( 0 ) differential cross section, mb per steradian 
0 scattering angle, radians 

X eikonal phase shift, dimensionless 

Subscripts ; 

C Coulomb 

cm center of mass 

c central 

s spin-orbit 

T target 
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PHENOMENOLOGICAL OPTICAL POTENTIAL ANALYSIS OF 
PROTON-CARBON ELASTIC SCATTERING AT 200 MeV 


Hari B. Bidasarla 
Old Dominion University 
Norfolk, Virginia 

and 

Lawrence W. Townsend 
Langley Research Center 
Hampton, Virginia 


SUMMARY 

Differential cross sections for 200 MeV protons elastically scattered 
from ^^C are analyzed utilizing a local, complex, spin-dependent optical 
potential with a harmonic v/ell radial dependence. Analyses were performed 
using the WKB and eikonal approximations. For the latter, first-order 
corrections to the phase shifts were incorporated to account for the 
spin-orbit contribution. Large disagreement between theory and experiment 
v/as observed when the usual Thomas form for the spin-orbit potential was 
utilized. Substantial improvement was obtained by allowing the parameters 
in the central and spin-orbit potential terms to vary independently- 
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INTRODUCTION 

The development of an accurate heavy ion (HZE) transport theory Is 
required to support cosmic ray shielding analyses for future space appli- 
cations. Transport theory development, in turn, requires a simple, yet 
accurate, HZE reaction theory as input. For the past several years, such a 
theory has been under development for incident projectile energies above 
several hundred MeV/amu (refs. 1-7). At present, no satisfactory theory 
exists for energies below (approximately) 200 MeV/amu where the elkonal 
(high energy) approximation appears to become inadequate (refs. 3 and 6). 
At lower Incident kinetic energies, contributions from low order partial 
waves (orbital angular momentum states) increase appreciably. In addition, 
the incident ion's trajectory deviates significantly from its initial 
straight line path such that the eikonal approximation may no longer yield 
accurate predictions. At these lower energies, a better approximation is 
the Wentzel-Kramers-Brillouin (WKB) method (ref. 8). This work describes 
our initial investigation into utilizing this approximation method to 
analyze elastic scattering. Proton-carbon scattering was chosen for 
analytic convenience due to its simplicity (elementary rather than ccwnpos- 
ite projectile) and because it is the subject of much current theoretical 
and experimental interest (refs. 9-11). Additionally, a harmonic well 
potential parameterization rather than the usual Woods-Saxon form (refs. 9 
and 10) was chosen since it more accurately represents the actual carbon 
nuclear distribution (refs. 7 and 12) and should therefore yield agreement 
with experimental data for fewer free parameters than the 24 required by 
the double-tJoods-Saxon potential analysis of reference 9. 
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ANALYTIC METHODS 

Optical Potential Parameterization 

The optical potential chosen for this analysis consisted of a central, 
Vc, and spin-orbit, Vg, contribution (ref. 8) 

V (r) » (r) + S • t Vj (r) (1) 

where o Is twice the spin angular momentum, s, and L the orbitaj 
angular riKmentum vector for the Incident proton. The scalar product o • L 
Is 


o 


t = 


L 


, parallel 


-a + 1) 


antiparallel 


( 2 ) 


where parallel/antiparallel refers to the relative orientation of a 
with respect to the total angular momentum vector 


^ = t + (1/2) a . (3) 

From reference 13, the spin-independent central potential for p - 
scattering, using a harmonic well charge distribution for ^,s of the 
form 


► 


Vj (r) ■ Aq (1 + AjP^) exp (-A^r^) (4) 

where Aq is complex and Ai and A 2 are real. The usual Thomas- 
Ferml form for the spin-orbit contribution is (ref. 11) 

V- (r) = 1 d V_ (r) (5) 

^ r ^ ^ 

which can also be parameterized as 

Vj <^) = (1 + Bjr^) exp (-B^r^) (6) 

with Bq complex and Bi and B 2 real. Using equations (4) and (5) It 
can be shown that the Bi (i = 0 -2) are functions of the Ai and 

that B 2 = A 2 . In this work, tnese parameters are treated as completely 
free. Hence, we find that the restrictions Imposed on the Bi by 

equation (5) do not hold and in particular, B 2 ^ A 2 (see RESULTS). 
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Scattering Formalism 

From reference 8, the scattering amplitude is 

F (0) » f (0) + o • n g (0) (7) 

where n is the unit vector normal to the scattering plane and 

f (0) = f (0) + (2ik)“^ i exp (2in, ) (8) 

^ L=0 

X |(L + 1) [exp (2i6j^"^) - 1] 

+ L [exp (2i5^^“) - 1] | P|_ 

00 

g (0) = (2k)"^ I [exp (2i<5 - exp (2i6 ")] Pj^ (cos 0) (9) 

L=1 *- 


where fc(0)» the coulomb scattering amplitude is 

(0) = -[6/2k sin^(0/2)] exp [2in^ - i8 In sin^(0/2)] . (10) 

In equations (8) through (10) above, k is the center of mass wave number, 
the Si are the phase shifts (to be determined in the following 
sections), the Pl(cos 0) are Legendre polynomials, and Pl^(cos 0) 
are associated Legendre functions of the first kind (see ref. 14). 
Additionally, for the coulomb contribution, we use 

8 = Z^e^/hv (11) 

Hq = Arg r (1 + i8) (12) 

and 

=* + tan"^ (8/L) (13) 

where v is the relativistic relative velocity between the projectile and 
target. 

The elastic differential cross section is then given by 

do - 0 (0) = I F (0) (14) 
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which, frcn equctiCf i (7) through (10) yields 


a 

where 


(0) - I t (0) 1^ I g (0) f * 2 Re (f*g) <n • 
Pi Is the polar izatlon of the Incident beam. 
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Elkonal Phase Shifts 

In terms of the elkonal phases, x» the partial wave phase 
first order corrections (rof. 11) to account for the 
contribution, are 

5|_* • (1/2) Ixj (5) ± (kh ? 1/2) X, (S)} 


where 


kb = L + 1/2 


and 


(S) = -(m/k) / (r) dz 

mCO 

Xj (B) = -(m/k) / (r) dz . 


shifts, with 
spin-orbit 

(16) 

( 17 ) / 

(18) 

(19) 


I 
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WKB Phase Shifts 

To first order, the WKB phase shifts for outgoing scattered waves are 
given by (ref. 15). 


«f • / k [1 - U*(r,L)]^^^ dr 


- / k [1 - (r.L)]!'^ dr 


where 


L(L + 1) ^ 2,e' 

Ur (r,L) p- 

^ 2mr^E Er 


and 


IT (r,L) - Vg^^/E 


( 20 ; 


( 21 ) 


( 22 ) 


with 



L(L + 1) 



\ * '■S 


Vc - a * 1) V, 



(23) 


In the above expressions, E is the projectile kinetic energy In the 
center of mass system. The complex turning points r|_* and r^ are 
those values of r for which 


1 - U~ (r,L) = 0 


and 


1 - (r.L) = 0 


(24) 

(25) 


respectively. These were located using a computer search routine based 
upon the Newton-Raphson method in the complex plane (ref. 16). The Inte- 
grals for the complex phase shifts in equation (20) were then evaluated 
numerical ly. 
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RESULTS 

Figure 1 displays theoretical cross sections obtained using 
equation (15) with elkonal phase shifts obtained from equation (16). For 
cwnparlson, the experimental rosiilts frcrni reference 9 are also displayed. 
The polarization of the incident beam was approximately 70 percent 
(ref. 9). The optical potential parameter values utilized for the anal- 
ysis, listed In Table I, v/ere obtained using a gradient-search least- 
squares fit to minimize the chi-squared deviation between the theoretical 
predictions and experimental results (ref. 17). As displayed In Figure 1, 
the agreement betn^en theory and experiment is excellent except for small 
angles (©cm ^ 10*). This disagreement is not unexpected when using an 
elkonal formalism (ref. 18). Table II displays the calculated parmeters 
for the Thomas-Fermi spin-orbit form from equation (5). Canparing Tables I 
and II shows the calculated parancter values to be quite different than the 
pareaneter values obtained with the computer search. 

Figures 2 and 3 display results utilizing the VflCB approximation. 
Complex turning points obtained from equation (24) are shown In Figure 2. 
For these calculations the potential parameters used In equations (4) and 
(5) are those listed in Table I. Mote that all turning points lie on or 
near the real axis. Figure 3 displays the cross sections obtained using 
'•/KB phase shifts. Improved agreement between theory and experiment is 
noted for small angles. There Is, however, significant disagreement for 
large angles (0rm >55*). This is likely due to the use of the potential 
parameters obtained from the eikonal phase shift search. Improvement 
should be obtained when the parameter search using WKB phase shifts is con- 
ducted, since the previously mentioned eikonal shortcoming for small angles 
should be corrected. A parameter search utilizing WKB phase shifts 
will, however, be very tim-^-ennsuming since any change in a parameter 
requires that all turning points and phas shifts be re-determined. Upon 
completing the paraneto’" sca-'ch with the WKt-. method, the conclusion con- 
cerning the use of the Thrmas-'^e’-ni -‘or the spin-orbit rontribution 
should be rechecked since the previous conclusion may no longer be valid. 
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CONCLUDING REMARKS 

This work has demonstrated that reasonably good agreement between 
theory :>nd experiment c«n be obtained with only 8 parameters In the nuclear 
potential. The use of potential parameters, obtained from an elkonal 
formalism, for a WKB analysis. Is probably not adequte but may serve as a 
useful starting point for a WKB parameter search. Other Improvements to 
thA theory could be made by discarding the point-charge Coulomb potential 
111 favor of a more real1«t1i. extended charge distribution. Cou1oaA> spin- 
orbit effects and Pauli blocking effects should also be Included. The 
latter, however, will Increase the nimiber of free parameters from 8 to 16. 
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Table I. - Optical Potential Parameters from the Eikonal Formalism Search 


Aq = -4.95 - i 48.01 MeV 
Ai = 0.380 fm-^ 

A2 = 0.495 fm-^ 

Bo = 8.52 + i 1.81 MeV 
Bi = 0.431 fm-^n 
B2 = 0.688 fm-^ 
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Table II. - Calculated Thomas-Fenni Spin-Orbit Parameters 

Bo = 1.14 i 11.04 MeV 
Bi * 1.636 fm-2 
02 = 0.495 fm-^ 
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Figure 1. 200 MeV proton -carbon elastic scattering 
differential cross sections. The theoretical curve 
was obtained using eikonal phase shifts. 
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Figure 3. 200 IVteV proton-carbon elastic scattering 
differential cross sections. The theoretical curve 
was obtained using WK6 phase shifts. 




